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Abstract 

We discuss the application of the resolvent technique to prove stabil- 
ity of plane Couette flow. Using this technique, we derive a threshold 
amplitude for perturbations that can lead to turbulence in terms of the 
Reynolds number. Our main objective is to show exactly how much con- 
trol one should have over the perturbation to assure stability via this 
technique. 
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1 Introduction 

We discuss stability of plane Couette flow via the resolvent method. Applying 
this method, one can derive lower bounds for the norms of perturbations that 
can lead to turbulence. Our aim is to discuss and clarify a point that has been 
overseen so far, which is to determine how much control over the perturbations 
one should assume to derive the stability result via the resolvent technique. We 
discuss the two dimensional case, but all the considerations can be used for 
the three dimensional case with minor technical changes. The main difference 
between the two and three spatial dimensions is that different resolvent estimates 
hold for each case, leading to different thresholds. This point will be made clear 
later on. We begin describing the problem and discussing previous works using 
the resolvent method. 

2 The problem and known results 

We are interested in the following initial boundary value problem: 

ut + (u ■ V)u + Vp = —Am 
K 

V • u = 

u(x,0,t) = (0,0) (i) 
u{x,l,t) = (1,0) 
u(x,y,t) = u(x + l,y,t) 
u(x, V, 0) = f(x,y) 
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where u : R x [0, 1] x [0, oo) — ► K 2 is the unknown function u(x,y,t) — 
(ui(x,y,t),U2(x,y,t)). The positive parameter R is the Reynolds number. The 
initial condition f{x,y) is assumed to be smooth, divergence free and compati- 
ble with the boundary conditions. The pressure p{x, y, t) can be determined in 
terms of u by the elliptic problem 

{Ap = -V • {{u ■ V)w) 
p y (x, 0, t) = —u 2y y(x, 0, t) 
p y (x, l,t) = -^u 2! yy(a;, 1,*)- 

It can be easily seen that [/(a;, y) = (y, 0), P = constant is a steady solution of 
problem Q). The vector field U(x,y) = (y,0) is known as Couette flow. 

Using the resolvent technique, one can prove and quantify asymptotic sta- 
bility for this flow. By quantification we mean the derivation a number M{R) 
such that disturbances of the flow with norm less than M(R) will tend to zero 
as time t tends to infinity. In other words, deriving a lower bound for the norm 
of perturbations that can lead to turbulence. For general discussion about the 
resolvent technique, see [H] and 0]. 

This problem has been studied for the 3 spatial dimensions case in [5 , and a 
threshold amplitude for perturbations was found to be of order 0(R~~). The 
estimates of the resolvent of the linearized equations governing perturbations 
where those found numerically in [8] and JU| , predicting the resolvent constant 
of the linear operator associated with the problem to be proportional to R 2 . 
In |3, the resolvent technique was used again to prove the stability of the 3 
dimensional problem but the estimates for the resolvent constant were those 
in [2| . By using modified norms, the authors achieve M(R) of order 0(i? -3 ) 
for two of the components of the perturbation, and of order 0(R~ A ) for the 
remaining component. Our approach uses again the resolvent technique, and 
we use the same norms as [H], with the obvious modifications for the 2 spatial 
dimensions case. We show that this approach leads, in our case, to a threshold 
amplitude of order 0(R~ 3 ). We note that our argument is the same used in 

with some minor differences. The only reason for the better exponent in 
our case is the better dependence of the resolvent constant on R for the 2 
dimensional case. In this case, the resolvent constant is proportional to R, as 
found in pQ. We carry out the argument in details again only because it is 
important for our aim, which is to clarify a subtle point that has been overseen 
in previous works: In [S], it was said that one needs control over the sobolev 
norm H 2 of the perturbation to assure stability. Later on, in [H], the authors 
note that the H 2 is not enough, and claim that one needs control over the 
norm H 4 . Actually, this is not enough yet, since in one of the directions, one 
needs control over six derivatives of the perturbation. This necessity is due to 
the pressure terms appearing in the problem. In section [S] we show in details 
estimates for these terms, and clarify the reason of this requirement. Moreover, 
our argument shows that derivatives of different orders of the perturbation scale 
differently with the Reynolds number. In other words, to assure decay of the 
perturbations via the resolvent method, one should require the perturbation to 
be small in some weighted norm involving six derivatives, where the weights 
depend on the Reynolds number R. 
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This work is divided in 4 sections: in section [3] we introduce some basic 
notation and derive the equations for perturbations of the Couette flow; in 
section 0] we derive estimates for the solution of the linearized equations for the 
perturbations; in [31 we use those estimates to prove asymptotic stability for 
the flow, and to derive the threshold amplitude M(R). In section [HI we derive 
carefully the estimates for the pressure terms involved in the problem. 



3 Notation and equations for the perturbations 

We denote by (•, •} and |j • || the L2 inner product and norm over 17 = [0, 1] x [0, 1]: 

{u,w) = / u-wdxdy ; ||m|| 2 = (u, u). 
Jn 

All the matrix norms that appear in this paper are the usual Frobenius norms. 
The usual sobolev norm of u over Q is denoted by 

n 

IMk»(n)=Ell^«ll 2 
3=0 

where D J denotes the j-th derivative of u with respect to the space variables. 
Unless stated otherwise, all norms in the space variables will be calculated over 
fl and therefore we will write || • || jjnrn) as || • || H n ■ We make use of a 2 dimensional 
version of the weighted norm || • ||^ used in Kreiss 5 : 

hll|=hl| 2 + ^II^H 2 + ^ll^l| 2 - (3) 
We also define another weighted norm || ■ \\ H e by 

IMIIre, = lkll^ + ^ll^ll 2 + ^P 4 "ll 2 

+ ^4 \\ u 2xxyyy || + \\ u 2yyyyy \\ + \\ u 2yyyyyy \\ ; (4) 

where u = (u\, u 2 ). 

The maximum norm over Q is denoted by | • loo. The norm || • ||^ is related 
with the maximum norm by the sobolev type inequality (see [2], Appendix 3, 
Theorem A.3.14) 

\-L<CR\\-\\%. 

Since we are interested in functions which are also dependent on time, we use 
that 

\u(;t)\l<CR\\ U (;ml, Vt>0, (5) 

where C is a constant independent of any of the parameters. 

We are interested first in proving asymptotic stability for the Couette flow, 
which is a stationary solution of that is, to prove that perturbations of 
the stationary solution that are small enough in some norm will tend to as 
t tends to infinity. More specifically, we will show that perturbations having 
norm || • \\ H e of order R~ 3 decay with time. 
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To this end, let U = U(x,y), P = P{x,y) be a stationary solution of 0. 
We can obviously use the Couette flow, but we think that the structure of the 
argument is easier to be understood if one uses any stationary solution. This 
will not change the estimates we will prove. We derive the equations satisfied 
by perturbations of this base flow. Let u(x,y,t) , p(x,y,t) be a solution of 
|QJ with initial condition f{x,y) = U(x,y) + ef'(x,y), where /' is divergence 
free and ||/'|| j^e ^ = 1. Then, e defines a unique perturbation amplitude. 
Write u(x,y,t) = U(x,y) + eu'(x,y,t) and p(x,y,t) = P(x, y) + epi(x, y, t) + 
e 2 p 2 (x 1 y,t). Then u',p'i,p'2 satisfy the system 

u' t + (u' ■ V)U + (U ■ V)u' + Vp[ + e(u' ■ V)u' + e\/p' 2 = -Au' 

R 

V • v! = 
u'(x,0,t) = (0,0) 
u'(x,l,t) = (0,0) 
u'(x,y,t) = u'(x + l,y,t) 
u'(x,y,0) = f'(x,y). 

The functions p\ and p' 2 are given in terms of u' by 

Ap[ = -V • ((u' ■ V)U) - V • ({U ■ V)u') 
1 



p' ly (x,0,t) 



R U 2yy 

I 



(x,0,t) 



and 



Ap' 2 = -V • ((«' ■ V)u') 
p' 2y (x,0,t) = 
p' 2y (x,l,t)=0. 

As we show in section \R\ the functions p[ and p' 2 can be estimated in terms of 
v! by 

||v P i(-,-,t)ll 2 < c(\\u\-rMm + ^h' 2yy (-rM + ^h' 2yy y(-r,t)\^ ,vt>o, 
l|vp' 2 (-,.,t)|| 2 < \\( u '-vy(;;t)\\\yt>o. 

From now on, to simplify the notation, we drop the ' in the equations above, 
and just write it, pi, p 2 . With this notation, the equations above are 



u t + (u- V)U + (U ■ V)it + Vpi + e(u ■ V)u + eVp 2 = 4 Au 

R 

V-u = 
u(x,0,t) = (0,0) 
u(x,l,t) = (0,0) 
u(x,y,t) = u(x+ l,y,t) 
u(x,y,0) = f(x,y), 

' Ap x = -V • ({u ■ V)U) - V • ({U ■ V)u) 
p ly (x,0,t) = —u 2yy (x,0,t) 

p ly (x,l,t) = —U 2 yy{x ) l,t) 



(6) 



(7) 
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and 

Ap 2 = -V • ((u • V)u) 

p 2y (x,0,t) = (8) 

P2y(x, l,t) = 0. 

Note that pi depends linearly on u. Moreover, for all t > 0, we have 
l|Vpi(-,-,t)|| 2 < C^( v ,t)||^+^|| U2TO (.,.,0|| + ^|| W2 ^(.,.,t)|^(9) 

||Vp 2 (.,-,i)|| 2 < \\(U-V)u(;;t)\\ 2 - (10) 

When the initial data is divergence free and the terms of pressure are given by 
the equations (JZJ) and (JHJ above, the solution u of problem © remains divergence 
free for all time t. Therefore, we drop the continuity equation and write problem 
© as 

u t = Cu — e(u ■ V)tt — eVp 2 
u(x,0,t) = (0,0) 

u(x,l,t) = (0,0) (11) 
u(x,y,t) = u(x + l,y,t) 
u{x,y,0) = f(x,y), 

where £ is a linear operator depending on the parameter R, defined by 

£u = -Att-(u-V)£/-(£/-V)u-Vpi, (12) 
R 

with pi given by Q . It is very important to note that this linear operator has 
also an integral part, which is the term Vpi. Moreover, as inequality @ shows, 
to estimate ||Vpi|| one needs three space derivatives of the second component 
of u, at least in one of the directions. 

This integral part of the operator C seems to be the point which was so far 
overseen. If one neglects that and consider only to the differential part of C, one 
will be led to conclude that || • ||#2 will be enough to estimate the right hand 
side of JI2J). 

We first apply the resolvent technique to prove stability of the stationary 
flow. For that end, it is convenient to have homogeneous initial conditions. 
Therefore, we transform the problem (|llfl to a similar problem with homoge- 
neous initial condition by defining 

v(x, y, t) := u(x, y, t) - e^f{x, y). (13) 

Note that v and u have the same behavior as t — > oo. Moreover, v given by l|13f) 
satisfies 

v t = Cv - e{{v ■ V)v + e~\v ■ V)/ + e"*(/ • V)w} - eVp 2 + F(x, y, t) 
v(x,0,t) = (0,0) 
v(x,l,t) = (0,0) 
v(x,y,t) = v(x + l,y,t) 
v(x,y,0) - (0,0), 

(14) 

where F(x,y,t) = e- i ((£+/)/-ee-*(/-V)/). Note that F, F t 6 L 2 ([0, oo); i 2 (0)), 
that is, both ||F(-, -, t)\\ 2 and \\F t (-, •, t)|| 2 are integrable over [0, oo). 
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4 Linear Problem 

We first consider the general linear problem 

v t = Cv + F(x,y, t) 
v(x,0,t) = (0,0) 
v(x,l,t) = (0,0) 
v{x,y,t) = v(x + l,y,t) 
v(x,y,0) = (0,0), 



(15) 



where ||T(-,t)|| 2 and \\F t {- 



t)\\ 2 integrable over the domain [0, oo): 
\F(;t)\\ 2 + \\Ft(;t)\\ 2 )dt<™. 



In our case of two spatial dimensions, resolvent estimates were found in pQ: 



R-,s)f < CiR 2 \\F(;s)\\ 



Res > 0, 



(16) 



where ~ stands for the Laplace transform with respect to t , s is its variable and 
C\ is an absolute constant, that is, it does not depend on any of the parameters 
or functions. One can prove, as in [5], Appendix A, that <|16f) implies 



\\v(.,s)\\%<CR 2 \\F(., S ) 



(17) 



where C depends on C\ and on U and its first derivative. Since for our problem 
U is fixed as the Couette flow, C is an absolute constant as well. From now 
on, we will use C for any absolute constant, and replace its value as necessary 
keeping the notation C . No attempt is made to optimize those constants, since 
the most important result is the dependence of the threshold amplitude on the 
Reynolds number. 

Using Parseval's relation, inequality l|17Jl for the transformed functions is 
translated to the original functions as 



\\v(;t)\\%dt<CR 2 



\F(- 



l dt. 



(18) 



pT poo 

Obviously, / \\v(-, t)\\ 2 ~dt < / \\v(-, t)\\ 2 ~dt , VT > 0. Moreover, since the 

Jo Jo 
solution of the equation up to time T does not depend on the forcing F{x,y,t) 

for t > T, we have 



K 



\\dt < CR 2 



\\F(;t)fdt 



VT > 0. 



(19) 



For our argument, we also need similar estimates for Vt- To this end, differentiate 
equation (fH)|l to get 



v tt = Cv t + F t (x,y,t) 
v t (x,0,t) = (0,0) 
v t (x,l,t) = (0,0) 
v t (x,y,t) = v t (x + l,y,t) 
, vt(x,y,0) = F(x,y,0) =: g(x,y), 



(20) 



G 



I) 



that is, v t satisfies an equation of the same type as ifPo)) . but with non-homogeneus 
initial conditions g(x, y) — F(x, y, 0). Performing the same type of initialization 
as before, that is, defining ip :— vt — e~ l g, we get a similar problem for ip, with 
homogeneus initial conditions and an extra forcing term. Using the estimates 
for the resolvent, and writing those in terms of v t , we get 

T \\v t (;t)\\%dt < \\F(x,y,0)\\% + CR 2 \\(C + I)F(x,y,0)f 

,T (21) 

+ CR 2 / \\F t {-,t)\\ 2 dt , VT>0. 
Jo 

Combining H19|) and l|21[l gives, for v the solution of 115|) . 

^ '\\v(.,t)\\% + \\v t (-,t)\\l)dt< \\F(x,y, 0)||| +Ci; 1 \l£+ I)F(.,:,,A))\\- 

+ CR 2 f (||^(-,t)|| 2 + ||F t (.,t)|| 2 )^ , VT>0. 
Jo 

(22) 

Now, using these estimates for the solution of the linear problem, we can prove 
a stability result for the nonlinear equation. 

5 Stability for the Nonlinear Problem 

The nonlinear problem 114|) is 

v t = Cv - e{(v ■ V)v + e- f {v ■ V)f + e"*(/ • V)v} - eVp 2 + F(x, y, t) 
v(x,0,t) = (0,0) 
v(x,l,t) = (0,0) 
v(x,y,t) = v(x + l,y,t) 
{ v(x,y,0) = {0,0), 

(23) 

where F(x,y,t) = e *((£ + /)/ — ee '(/ ■ V)/). We prove the following: 

Theorem 5.1 There exists eo > 7 eo = ea(R), such that if < |e| < eo, then 
the solution v{x,y,i) of \2ty satisfies 

lim |z;(-,t)|oo = 0. 

t — >oo 

Moreover, e = 0(Rr 3 ). 

Proof: We consider problem (|23|l as a linear problem with forcing 

G(x, y, t) := F(x, y, t) - e{(v ■ V)w + e~*{v ■ V)f + e -*(f ■ V)v} - eVp 2 . (24) 
Applying inequality (|22|l with forcing term G gives 

T (\\v(-M\+\\<-M\)dt< \\G(x,y : 0)\\l + CR 2 \\(C + I)G(x,y 1 0)\\ 2 + 



T 



CR 2 (\\G(;t)\\ 2 + \\G t (;t)\\ 2 )dt VT>0.(25) 







7 



From the definition of G, we have 



(\\v(-Ml + \\M-Ml)dt<2\\F(x,y,0)\\l + 2e"\\yp 2 ( Xl y,0)\\l 
+CR 2 \\(Cr + l)F{x, y, 0) || 2 + CR 2 \\ (C R + l)p 2 (x, y, 0) || 2 
+ CR 2 [ (\\F - e{(v ■ V)v + e _t (u • V)/ + e"*(/ • V)v} - e\7p 2 \\ 2 ) dt 

+ CR 2 [ (\\(F - e{(v • V)» + e-*(« • V)/ + e"^/ • V)v} - eVp 2 ) t \\ 2 ) dt. 
Jo 

(26) 

Since p 2 is given by JSJ, we have (see sectional 

||V P2 || < ||(«.V)«|| ; [[(Vi*)t|[<||((u-V)u)t||. 
Thus, using we can estimate Vp 2 by / and v. Moreover, 

l|Vp 2 (.,.,o)|| 2 <||( w .vM.,.,o)|| 2 ^||(/.v)/|| 2 , 

and since ||/||^ 6 = 1, inequality (|2*H|l gives 

i-T 



CR 2 / (||^|| 2 + ||F t || 2 )di + Gi? 2 e 2 / (||(».VH 2 + ||( Ut -VH| 2 + ||(7;-V) Ut || 2 )dt 
Jo T Jo 

CR 2 e 2 f (||e-*(« • V)/|| 2 + \\e~\f • V)«|| 2 + ||e-*(«* • V)/|| 2 + || e -*(/ • VM| 2 ) dt. 



Since 

F(a;, y,t) = e -*((£ + 7)/ - ee~\f • V)/), 

we have F(x,y,0) = (C + I)f - e(/ • V)/ := P/. With this notation, the 
inequality above is 



(\\v(;t)\\l + ht(;t)\\%)dt< \\rf\\\ + CR 2 \\{C + I)Vf\\ 2 + CR 2 \\{C + I)f\\ 2 
+ CR 2 e 2 \\(f ■ V)/|| 2 + CR 2 e 2 f (\\(v ■ V)«|| 2 + \\(v t ■ V)v\\ 2 + \\(v ■ V)v t \\ 2 ) dt 

T J ° 

+ CR 2 e 2 [ (||e-*(« • V)/|| 2 + [[ e -*(/ • V>|| 2 + \\e^(v t ■ V)/|| 2 + +\\ e - t (f ■ X7)v t \\ 2 ) dt. 



. (27) 

To apply the resolvent method, one needs control over the terms depending on 
/ of the right hand side of inequality (|27[) . Its second term is 

HOC + 7)P/|| = HOC + /)((£ + /)/ - e(/ • V)/)||. 

To estimate ||£ 2 /|| it is necessary to have, at least in the y direction, control 
over six derivatives of f 2 , the second component of /. The reason is that three 
derivatives of f 2 in the y direction are necessary to control Cf, due to the 
integral part of the operator C. Moreover, derivatives of different orders of the 
perturbation may have different scales with respect to the Reynolds number. 
These facts are shown in details in section |SJ As already mentioned, clarifying 
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and showing these estimates in details is our aim, since it is a point that has 
been overseen in previous works, and lead to mistakes about the necessary 
assumptions on /. 

We now continue the proof of stability. As mentioned above, since ||/||^6 — 
1, we can replace all the terms depending on / by an absolute constant and 
write inequality (|2*7|) as 

£ (|MI| + IHI|) dt < CR 2 + CR 2 e 2 jf \\(v ■ V)v\\ 2 dt 

+ CR 2 e 2 f (\\(v t ■ V)v\\ 2 + \\(v VMI 2 ) dt (28) 
Jo 

+ CR 2 e 2 f (We^iv ■ V)/|| 2 + || e -*(/ • V>|| 2 ) dt 
Jo 



+ CR 2 e 2 / (\\e-\v t ■ V)/|| 2 + || e -*(/ • VH| 2 ) dt. 
Jo 

From now on, we fix the constant C. To finish the proof, we use the following 
Lemma, which is proved later: 

Lemma 5.1 There exists eg > 0, eo = 0(i?~ 3 ), such that if < e < €q then 

£(\\v(-,t)\\ 2 ~+\\v t (- : t)\\l)dt<2CR 2 ,VT>0. (29) 
Now, using (JSJ and a simple one dimensional sobolcv inequality, we have 
max K-,^ <CRm^Jv(-,t)\\l<CR^l + ^—^j J ^ \\v(;t)\\%dt+ J ^ \\v t (;t)\\%dt. 
This implies 

S up\v(;t)\l <CR f~ (\\v(-M% + \\M-M%)dt. (30) 



a<t 

Note that in view of Lemma l5~Tl the right hand side of inequality l|30|) is finite. 
Letting a — > oo in (|30|l . we have that lim ^(-(t)^ = 0, which proves the 

t — >oc 

theorem. 

Proof of Lemma 15.11 First, note that for T > small enough, inequality 
(|29|l obviously holds. Now, suppose it does not hold for all T > 0, that is, there 
exists T > such that 



t 



(lK- J *)||| + IM-,i)ll|)^ = 2Ci2 2 . (31) 
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Using (J2SJ, we have: 

2CR 2 = £° (\\v\\% + \\v t f s )dt<CR 2 
+ CR 2 e 2 [ T ° (\\(v ■ VH| 2 + ■ V)vf + \\(v ■ X7)v t \\ 2 ) dt 

To ° 

+ / ° (\\e~\v ■ V)/|| 2 + ||e-*(/ • VH| 2 + \\e-\vt ■ V)/|| 2 + \\e~\f ■ V)v t \\ 2 ) dt. 
Jo 

(32) 

We now estimate the integrands on the right hand side of inequality (|32() by the 
integral on its left hand side. To this end, we will use the inequalities (|5jl and 

max |K,*)II|< f T " (\H;t)\\l + \\v t (;t)\\%)dt = 2CR 2 . (33) 

Since H|| = \\v\\ 2 +^\\Dv\\ 2 +^\\v xy \\ 2 , we have \\Dv\\ 2 < R\\v\\ 2 ~. Therefore, 
for each < t < T : 

\\{(vV)v}(;t)\\ 2 < H;t)\UDv(;t)\\ 2 (34) 

< (CR\\v(;t)\\%) (R\H;t)\\%) < 2CCR 4 \\v(;t)\\% 

||{(«fVM(.,t)|| a < \v t (;t)\UDv(;t)\\ 2 (35) 

< (CR\\v t (;t)\\%) (RM;t)\\%) < 2CCR i \\v t {;t)\\% 

\\{(vV)v t }(;t)\\ 2 < Ivi-rftWDvti-M 2 (36) 

< 



CR\\v(;t)\\s) {R\\vt{;t)\\%) <2CCR i \\v t {;t)\\%, 

lie" 4 {(« • V)/} (-,t)|| 2 < e- 2t \v{;t)\lo\\Df\\ 2 < CR\\v(;t)\\l, (37) 

|| e - 4 {(/ • V)v} (;t)\\ 2 < e-*\f\UDv(;t)\\ a < R\H;t)\\% (38) 

lie" 4 {(vt ■ V)/} (-,t)\\ 2 < e- 2t \v t (;t)\l\\Dff < CR\\v t {-,t)\\% (39) 

lie"* {(/ • VK} M)|| 2 < e- 2t \f\UDv t (-,t)\\ 2 < (40) 
Applying ©, E3, dSHJ, <E1, © to 63 gives 

2Ci? 2 < Ci? 2 +Ci? 2 e 2 jecoR 4 ^ ° (|K-,t)||!+ ||«t(-, *)[[!) dtj ^ 
= CR 2 + CR 2 e 2 {l2C*C 2 i? 6 } . 
This implies 

1 < 12CC 2 i? 6 e 2 , (42) 



10 



which is equivalent to 



e > 



1 



where K := cV 12C. Therefore, if e < 
proves the Lemma. 



CR 3 Vl2C KR3 
1 



(43) 



, equality lj3T|) never holds. This 



6 Estimates for the pressure terms 

We now turn our attention to our main objective, which is to show why one 
needs to assume control over at least six derivatives of f 2 in the y direction to 
apply the method above. This necessity follows from estimates (0 and l|lUfl 
for the pressure terms Px(x, y, t) and p2(x, y, t). So, we begin by showing these 
estimates. 

Theorem 6.1 If pi(x,y,t), p2(x,y,t) are the solutions of 



and 

then 
||Vpi( v ,i 

||Vj>2(v,l 



Api = -V • ((u ■ V)U) - V • ({U ■ V)u) 

P ly (x,0,t) = —U 2 yy(x,0,t) 
p ly (x,l,t) = -=U2yy(x,l,t) 

K 



Ap 2 = -V • ((it • V)m) 

P 2 y(x,0,t) = 
P2y(x,l,t) = 0, 



(44) 



(45) 



< 
< 



Cn\u(;;: 

||(u-V)u(. 



1 



1 



Iff 1 + ^2lK OT (-,-,£)|! + -^2ll"2y TO (-,-,t)|[^6) 
*M|2 



(47) 



for all t > 0, where C is an absolute constant. 



Note that the inequalities above are for norms with respect to the space variables 
x and y. Therefore, to simplify the notation, we prove them for functions 
depending only on these variables. As before, — [0, 1] x [0, 1], and all norms 
are over J7. For clarity of the presentation, we separate the proof of Theorem 
Ift.ll into two Lemmas. 



Lemma 6.1 Let g 
function satisfying 



x [0, 1] 



2 , g( x ,y) = (gi(x,y),g2(x,y)), be a C° 



Ifh 



[0,1] 



g(x,l) =g(x,0) = (0,0) 
9{x,y) = g(x + l,y)Vx G 

is the solution of 

Ah = V • g 

h y (x,0) = 

h y (x, 1) = 
> h(x,y) = h(x + l,y), 



(48) 
(49) 



(50) 
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then 

liv/^ii 2 < || 5 || 2 . 

Proof: If A/i = V • g, then 

(h xx + h yy )h dxdy = (g lx + g 2y )hdxdy. 
i Jn 

Through integration by parts, 



(51) 



(h\ + hi) dxdy + 



- / g\h x dxdy 



h x hv x dS - 



on 



on 



h y hv v dS 



gihy dxdy 



g\hv x dS 



on 



g 2 hv v dS, 



on 



< 5 1|* 



| 2 + ^ll^|| 2 



1.92 I 



;\\h y \ 



where v x and v v denote the components of the outer normal to dVL in the x 
and y directions respectively. From the conditions satisfied by h and g at the 
boundary, the boundary integrals above vanish. Then, 

I'^ll 2 = / (h x + hy)dxdy= / g\h x dxdy+ / g 2 h y dxdy. 
Jn Jn Jn 

Using the Cauchy-Schwarz inequality, 

||v*|| a < IMIIIM + MIW 

This implies the desired estimate 

[ivftf < Nl 2 . 

The Lemma above gives inequality (|47|l for p 2 (x,y,t), the solution of l|45|l. 

As mentioned before, the estimates to be proved do not depend on the 
variable t. Therefore, we write simply u(x,y) = (ui(x,y),u 2 (x,y)), for u the 
solution of problem ©. We remind the reader that U(x,y) = (y, 0). We prove 
the following lemma, completing the proof of Theorem 16. II 



' ' |2 1 l\Nh\\ 2 . 



Lemma 6.2 Ifh:Rx [0, 1] 



is the solution of 



( Ah = -V • (u ■ V)U - V • (U ■ V)it 

K( x i°) = -^ u 2yy(x,0) 
hy{x,l) = -jjU 2yy {x,l) 

> h(x,y) = h(x + l,y), 



(52) 



then 



Vh\\ z < C [ \\u\\ J H1 + jp\\u 2yy \\ 2 + -^\\u 2yyy 



R 2 



(53) 
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(54) 



(55) 



Proof: We begin by noting that if hi, h 2 are the solutions of 

Ah! = -V ■ (u ■ V)U - V • (U ■ V)u 

h ly (x,0) = 

hi y (x, 1) = 

hi(x,y) = hi{x + l,y), 

and 

r a/i 2 = 

h 2 y(x,0) = —U2yy(x,0) 

H 

h 2 y(x,l) = ftU 2 y y (x,l) 

> h 2 {x,y) = h 2 (x + l,y), 

then ft, = + h 2 is the solution of problem {52} . Therefore, to prove 1)53(1 . it is 
sufficient to prove estimates for fti and /12, solutions of l|54l) and {55} respectively. 
For hi, lemma RTT1 implies 

||V/ii|| 2 < \\(u-V)U+(U ■ V)u|| 2 < 2|M| 2 + 2||L>w|| 2 = 2\\u\\ 2 H1 . (56) 

To prove the estimates for h 2 , expand in a Fourier series in the x direction. The 
Fourier coefficients h 2 (k,y) satisfy 

-k 2 h 2 + h' 2 ' = (57) 

h 2y (k,0) = -u 2yy (k,Q) (58) 
ix 

h 2y (k,l) = ^-u 2yy (k,l), (59) 
H 

where ' denotes the derivative with respect to y. Consider first the case k =/= 0. 
To simplify notation, let 

otk ■■= -j^u 2 yy{k,Q) 



13k ■= ^U 2 yy{k,l). 



Using a one-dimensional sobolev type inequality, we estimate a k and (3 k by 

1 C 
\a k \ 2 < ^max \u 2yy (k,y)\ 2 < (\\u 2yy (k, -)\\ 2 + \\u 2yyy (k, -)|| 2 ) (60) 

\Pk\ 2 < ^maxj^^y)! 2 < ^ (\\u2 yy (k, -)\\ 2 + \\u 2yyy (k, -)\\ 2 ) (61) 

where C is an absolute constant. As before, we keep the notation simple by 
using C to represent any absolute constant, whose value can possibly change for 
different inequalities. 

The general solution of the differential equation (|57H is 

h 2 (k,y) ^ a k e^~ 1 '> +b k e-^. (62) 



13 



Imposing the boundary conditions H58|) and we determine the coefficients 
a k and b k : 

e' fc l k 1 a k b 1 ft e' fc l Qfc 

e |fe| _ e -|fc| |^,| e |fe| _ e I i e |fe| _ e -|fc| e |fe| _ e -|fe| 



Therefore, 



, k 2 k 2 
Using (EOJ), (ED, E3, E3>, JMJ>, we have 

fc 2 |MMH 2 < Ck 2 (\a k \ 2 \\e) k ^\\ 2 + M 2 ||e- |fe| l 2 ) < j^(\a k \ 2 + \[3 k \ 2 ) 

< C(K| 2 + |/3 fe | 2 )<^(||S 2ra (fc,.)l| 2 + ll^(fc,.)l| 2 ) (65) 

and 

fh 2y (k, -)H 2 < Ck 2 (\a k \ 2 \\e^y-V\\ 2 + |6 fc | 2 || e -'^|| 2 ) < ^(\a k \ 2 + \? k \ 2 ) 

< C(|a fc | 2 + |/3 fc | 2 )<-J(||S 2ra (fc,.)l| 2 + ll^(fc,-)l| 2 )- (66) 

For the k = mode, we solve equation (|57fl under boundary conditions (|58f) . 
<|59|) directly. Note that the divergence free condition satisfied by u assures that 
this problem is solvable, as expected. Differentiating the solution with respect 
to y, we get 

h 2y (0,y) = —u 2yy (0, 1). 

Then 

HMO, Oil 2 < ^ max \u2y V (0,y)\ 2 < ^ (\\u 2yy (0, -)|| 2 + ||2 2wy (0, -)l| 2 ) ■ 

(67) 

Using TO and j67 |) . 

Iiv^n 2 = £(fc 2 ||/; 2 (MI! 2 + IIMMII 2 ) 

< ^Edl^^OII' + ll^^,-)!! 2 ) (68) 

C 

= (ll^all 2 + \\u2yyy\\ 2 ) ■ 

Therefore, by (|56|l and l|()8l) . we have that h = hi + h 2 , solution of H52|) . satisfies 
||V/i|| 2 < C (jmi 2 ^ + -^h2 ra || 2 + ^II^H 2 ) • (69) 
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This finishes the proof of the Lemma. This Lemma applied to pi(x, y,t) com- 
pletes the proof of theorem 16. II 

In section [5J we used that the terms depending only on / on the right hand 
side of (|2*7Jl can be estimated by C||/||^ 6 . This fact is a consequence of Theorem 
16.11 In fact, it is easy to see that the term that requires highest derivatives of 
/ to be estimated is ||£ 2 /|| 2 . Therefore, we should bound this term in a sharp 
way, that is, estimate it using a norm for / that involves derivatives of lowest 
possible order. We show this norm to be ||/||.ffe^. 

Recall that 

Cf = iA/ - (/ • V)U - (U ■ V)/ - Vp, (70) 

where U — (y, 0), and p is the solution of 

f Ap = -V-((/-V)t/)-V-((t/-V)/) 

ft/ 0,1) = ^hyy(x,l)- 

Using Theorem 16. II we get that 

HAfll 2 < C Uff m + ^\\f xx \\ 2 + ^WfyyW 2 + ^\\f2yyy\\^j ■ (72) 

Therefore, 

\\C 2 f\\ 2 < C (\\£f\\ 2 Hi + ^\\Cf xx \\ 2 + ^WQyyW 2 + ±- 2 \\{Cf)2yyy\\ 2 ^ . (73) 

Straightforward computations show that 

\\£f\\m < c(\\f\\ 2 H2 + -^\\D 3 f\\ 2 + -^\\D 4 f\\ 2 
-^\\£f xx \\ 2 < ^(ll/ll^ + P 3 /l| 2 + ^P 4 /l| 2 + ^ll/ 2ra ^|| 2 

-^\\£fyy\\ 2 < ^(ll/ll^ + p 3 /l| 2 + ^p 4 /l| 2 + ^ll/2, rara || 2 
-j^\\{£f)2yyy\\ 2 < ( || jlyyy || 2 + || f2 X yyy \ | 2 ) 

c 

The inequalities above, together with l|73|) . imply 



R 4 (Il/2i/y2/xx|| + ll/2yyyyy|| + Wf2yyyyyy\\ 



\C 2 f\\ 2 < c(\\f\\ 2 H2 + ±\\D 3 f\\ 2 + ±\\D 4 f\\ 2 + j-Jf 2xxyvy \ 



+ ~j^\\f2yyyyy\\ 2 + -j^\\f2yyyyyy\\ 2S J — C\\J\\ 2 H t^. (74) 

Note that this estimate is sharp in view of theorem l6.1l and the arguments above: 
one needs at least six derivatives of fa to be able to estimate ||£ 2 /|| 2 . 
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7 conclusions 



As showed in section [5] perturbations with norm || • \\ H e of order R~ 3 decay 
with time. Note that this also shows that even though one needs control over 
derivatives of high order of the perturbation, there are different scales for those 
derivatives. In fact, if 

\\f\\m m = ll/ll 2 ff2 + ^P 3 /ll 2 + ^P 4 /ll 2 + ^4ll/ 2 ™l| 2 

+ ^4 WhyyyyyW 2 + ^4 \\hyyyyyy\\ 2 = 0(R 6 ), (75) 

then ||/||^2 is of order i? -3 , ||-D 3 /|| and ||-D 4 /|| are of order i?~ 2 , and \\f2xxyyyW-, 
WhyyyyyW and WhvvvvvvW are of order Rr l . 

In applications, control over this many derivatives is too restrictive. One 
possible way to avoid this requirement is to incorporate a smoothing property 
of the system to the argument. We hope to address this question in the future. 
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